The activated dynamics exponent i/j for glasses 
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The activated dynamics exponent which relates the correlation length in glasses to the alpha 
relaxation time is argued to equal d, the dimensionality of the system. This result is obtained via a 
detailed study of relaxation processes in small systems of disks confined in a square and extended 
to large systems by employing the mosaic picture of glasses. It is in excellent agreement with the 
simulational studies on hard spheres of Berthier and Witten. 
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As the temperature T of a fragile glass is reduced, 
the alpha relaxation time r a , the fundamental glass time 
scale, rapidly increases. In the old theory of Adam and 
Gibbs [1], and the more recent random first-order transi- 
tion (RFOT) theory of glasses [2] , [3] and its elaboration 
as a mosaic picture (for a review see [4]), this increase is 
attributed to cooperatively rearranging regions (CRR) of 
size £{T). This length scale increases as the temperature 
decreases and has actually been determined numerically 
from a study of "point-to-set" correlations [5, 6]. It is 
also expected to be related to the linear extent of the 
dynamical heterogeneities [4]. £ and r a are thought to 
be related by the activated dynamics equation 

r a -roexp[S(e/^)], (1) 

where B is a non-universal constant, a is a length of the 
order of a molecular diameter and To is a microscopic 
time scale. Much interest attaches to the value of the 
exponent ijj. In this paper we shall argue that 

^ = d, (2) 

where d is the dimensionality of the system. This result 
should hold both for hard spheres and disks, and might 
also extend to soft potentials. For hard spheres (d = 3) 
or hard disks (d = 2), the time scales and length scales 
are increasing functions of the packing fraction <fi rather 
than the temperature [7]. In general, one expects that 
the correlation length grows as 

c ~ (w/O 7 ' (3) 

where in the approach of Fullerton and Moore [8] (see 
also [9]), 



is the spin glass domain wall energy exponent, and for 
d = 3 is approximately 0.2 [10-12]. max in Eq. (3) is 
close to (or even identical to) the random close packing 
fraction (/> rcp ~ 0.635. Then, if ^ = d : we predict that 
the alpha relaxation time varies as 

r a ~ r exp [B(l - , (5) 



where the exponent 5 = i\)v = 2d/(d — 26) ~ 2.3 in three 
dimensions. Berthier and Witten [7] found that 5 ~ 2.2 
in their study of hard spheres, when they used a value of 
0max = 0.635 for their fitting. Note that Eq. (4) applies 
to correlations which exist in the supercooled liquid up to 
alpha relaxation times. It does not apply to behavior on 
time scales very much longer than the alpha relaxation 
time, when real supercooled liquids may phase separate 
(if they are mixtures of species) or simply crystallize. 

The basic idea of Adam and Gibbs [1] and the mo- 
saic picture [4, 13] is that regions of size less than £ are 
glassy and rarely relax. (In the simulations of Ref. [6] , for 
cavities of a size larger than £ the particles could escape 
from their initial configurations; but for cavities of size £ 
or smaller, the particles were often trapped near their ini- 
tial configurations when realistic, non-swap Monte-Carlo 
dynamics were used.) It is by a careful examination of 
just how the atoms within a region of linear extent ~ £ 
have to move in order to allow the atoms to escape from 
their initial configurations that we are able to argue that 
ip = d. On the other hand, RFOT with wetting predicts 
that i/j = d/2 and v = 2/d [2] so that 5 = 1, which is the 
Vogel-Fulcher-Tamman law [5]. If one uses the RFOT 
form v = 2/d, with our value of ip = d, then 5 = 2, 
which is also the value obtained by the very different ar- 
guments of Schweizer [14]. The simulations in Ref. [5] 
were claimed to give ijj = 1, a value which is hard to 
reconcile with the value of 5 of Berthier and Witten [7] . 

We shall illustrate the process of escape from the initial 
configuration for two simple systems, consisting of either 
two or five disks confined in a square. At high density, 
escape involves squeezing through a narrow neck in con- 
figuration space, which for these simple systems we can 
make explicit. We shall show that transition state theory 
provides a quantitative account of the slow relaxational 
processes in these small systems. The transition state 
theory is then generalized to large systems of hard disks 
or spheres, when the length scale of the dynamical het- 
erogeneities is of order £, and we derive Eq. (2). (We shall 
assume throughout that the dynamical and static length 
scales are proportional when they are both large [15].) 
Our generalization from the explicit calculations where 
the number of disks, AT, is small (2 or 5) to the case of 



2 




FIG. 1. Mutually inaccessible configurations of a system of 
two disks, illustrated for r = 0.27 L, which is greater than the 
critical value r — 0.25 L below which the disks can pass each 
other. 



large TV involves assumptions about the nature and va- 
lidity of the mosaic picture [4], but the good agreement 
for the value S with the results of Ref. [7] suggests that 
the assumptions made are actually correct. 

We begin with some definitions for our study of two 
and five disks confined in a square: L is defined as the 
length of one side of the square; r is the radius of the 
disks; and the packing fraction <p is defined as the fraction 
of the area of the square that is covered by the disks, 
(j) = Nnr 2 /L 2 . A configuration of the disks is defined by 
the Nd coordinates of the disk centers (for disks d = 2); 
every configuration belongs to a state, which is the set 
of configurations that can be reached from it without 
violating the no-overlap constraint. 

Two- disk system. The system of two disks confined to 
a square has been considered previously by Speedy [16]. 
For r < L/4 (or <j) < 0i g ~ 0.3927) the disks can pass each 
other, though this becomes more difficult as r — » L/4. 
For larger values of r, the configuration space (disregard- 
ing the identity of the disks) is broken into two states in 
which (for r -> L/(2+\/2), or (j) g ~ 0.5390) the disk 
centers lie on the same diagonal of the square (see Fig. 1). 
These limiting configurations are the inherent structures 
introduced in Ref. [17]. Speedy [16] has considered the 
thermodynamics of this system and, in particular, finds 
weakly non-analytic contributions to the thermodynamic 
quantities, such as the pressure. Bowles et al. [18] have 
studied the dynamics of a similar system of two disks con- 
fined to a channel, and find that the mean first-passage 
times of the disks past each other varies as (1 — 0/0i g ) _Q! , 
where a = 2 for inert ial motion of the disks. Our event- 
driven molecular dynamics results (Fig. 2) confirm the 
value of the exponent, a = 2, for the case of two disks in 
a square box. 

The origin of this behavior can be obtained from tran- 
sition state theory [18, 19]. In this well-studied approxi- 
mation, which works best when the transition rate over a 
barrier is small, the transition rate R between two states 
varies as 

R=1/t~ Z x jZ, (6) 

where Z* is the partition function evaluated at the top 
of the barrier; see Ref. [19] for a full description of the 
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FIG. 2. Molecular dynamics transition rates (1/r) between 
met ast able "glassy" states of the two-disk system, as the pack- 
ing fraction, 0, approaches the critical value 0i g ~ 0.3927 cor- 
responding to r = 0.25 L. The gradient of the straight-line fit 
is 1.94, which is close to the expected value 2. 



transition state formalism and the definition of Z*. In 
the case of two disks passing this means that instead 
of the full partition function integral over (x\,yi) and 
(^2,2/2), there is a constraint that 2/1=2/2 = V, so it 
is effectively a three-dimensional integral. The integral 
over y gives a trivial factor (L — 2r) and the remaining 
two integrals give a factor (1 — (/>/</>i g ) 2 m the limit <p — > 
0i g . Z itself is essentially just a constant: it has a very 
mild singularity, Z reg + C(l — </>i g /(/>) 5 / 2 , when the packing 
fraction <fi approaches <p\ g from above [16]. Transition 
state theory thus predicts a slope of 2 in Fig. 2 for the 
dependence of the relaxation rate on packing fraction. 

Configurations of the five- disk system. The configu- 
ration space of five disks confined to a square has been 
analyzed previously: Bowles and Speedy [20] have dis- 
cussed the thermodynamics and dynamics; Hinow [21] 
has studied the jammed states of this system; and Carls- 
son et al. [22] have given a detailed analysis of how the 
topology of the configuration space depends on r. We re- 
fer to Fig. 3 for configurations of the disks at two critical 
values of the radius. 

Below the liquid-glass critical point, i.e., for r < r cg ~ 
0.1863 L, the system is fluid (any pair of disks can ex- 
change position), but for slightly greater values of r the 
configuration space is fractured into two states: a "crys- 
talline" state in which one disk is surrounded by the 
four others, confined near the corners of the box; and 
an "amorphous" state in which all five disks lie close to 
the walls of the box and are unable to change their order. 
Above r = r gg ~ 0.1942 L, the state fractures further into 
four "frozen" amorphous states of the kind illustrated in 
Fig. 3 (b), in which one disk is confined near a corner of 
the box. Above r = r g ~ 0.1964 L, the system can exist 
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FIG. 3. (a) Representative configurations of the "crystalline" 
(left) and "amorphous" (right) states of a five-disk system, 
connected by a transition state (center), shown here with 
the largest radius, r c± 0.1863 L, for which the transition be- 
tween crystalline and amorphous states is possible; (b) two 
frozen amorphous states (left and right) and the transition 
state (center) connecting them, shown with the maximum ra- 
dius, r ~ 0.1942 L, for which a transition between amorphous 
states is possible. 



only in the crystalline state. 

It may be noticed that r cg ~ 0.1863 L differs signifi- 
cantly from the value 0.1871 stated by Carlsson et al. [22]. 
We have been unable to find a path between glass-like and 
crystal-like metastable states that passes via the config- 
uration proposed in their paper. We find, moreover, that 
their proposed state with r ~ 0.1871 L is not a stationary 
point of the softened potential energy function E intro- 
duced in Ref. [22]: instead, it is a minimum of \VE\ 2 at 
which VE 7^ 0. On the other hand, we can show that our 
own configuration at r ~ 0.1863 L lies on a path between 
crystal-like and glass-like states and also that this config- 
uration corresponds very precisely to an ordinary saddle 
point of E. Such a reaction path is illustrated by an 
animation provided in the supplement to this paper [23] . 

Dynamics of the five-disk system. As for the case of 
two disks, an event-driven molecular dynamics algorithm 
[24] was used to simulate the motion of the five-disk 
system and calculate the mean time of passage between 
metastable states. The initial velocities of the disks were 
chosen to be of fixed magnitude, but with random direc- 
tions. 

A very simple method was used in our work to iden- 
tify when a transition had taken place. For the transi- 
tion from crystal to amorphous states, the simulation is 
started in a typical "crystal" configuration with one disk 
[Fig. 3 (a), shaded] near to the center of the box. The 
shaded disk's first collision with any wall is an unambigu- 
ous sign that the transition to the amorphous state has 
occurred. Transitions between metastable amorphous 
states can be identified in a similar way, using the fact 
that every disk remains near the perimeter of the box. 
From Fig. 3 (b), we can see that a transition has occurred 




-2.8 -2.6 -2.4 -2.2 -2.0 -1.8 -1.6 

ln(l - 4>/4>cg) 



FIG. 4. Molecular dynamics transition times: (a) between two 
amorphous states; (b) between the crystalline and amorphous 
states. Both graphs have a noticeable curvature: the steepest 
gradients are approximately (a) —9 and (b) —6.5. 



if any disk makes a collision with a wall other than the 
one (or two) it was close to in the initial configuration. 
For each kind of transition, the time of first occurrence 
of the diagnostic event is recorded and the simulation 
restarted with random initial velocities. 

Results from our simulations are plotted in Fig. 4. 
They show that, as expected, the mean time of pas- 
sage increases rapidly as <ft approaches a critical value 

= ^cg o r ^gg- Because of the slowing-down of the dy- 
namics, our simulations were unable to reach a regime 
where a power-law dependence of r on (1 — (j)/4> c ) could 
be demonstrated unambiguously: the steepest gradients 
in the graphs of Fig. 4 (a) and (b) are approximately 
—9 and —6.5, respectively. The first of these results (for 
the transition between metastable amorphous states) is 
easy to understand in the transition-state theory. If the 
reaction coordinate is fixed at the value it takes in the 
transition state, the hard-disk constraints define a con- 
figuration space with nine spatial dimensions: accord- 
ingly, the constrained partition function (and hence r^ 1 ) 
should be expected to vary as (1 — (j)/(j) gg ) 9 ; we have veri- 
fied this dependence by a separate Monte Carlo sampling 
of the partition function. A similar argument could, of 
course, be made for the transition between metastable 
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crystal and amorphous states, if <p was sufficiently close 
to (j) cg ; our results from molecular dynamics suggest that 
the anticipated power-law behaviour is restricted to a 
very narrow range of <\> and that this asymptotic regime 
was not reached in our simulations. 

Our conclusion from these studies is that a transition 
between states in a region containing N particles gener- 
ally requires coordinated motion of all the N particles 
in order to squeeze through the neck in the phase space. 
The rate at which this will occur is, according to transi- 
tion state theory, 



1/t~(1/t )(1 



Neck)^ V(l 



Nd 



(7) 



where <\>j is the packing fraction of the jammed state clos- 
est to the initial configuration of the particles and </>Neck 
is the highest packing fraction below <pj at which a neck 
first opens to allow escape from the jammed state. In 
studies of jammed states in large systems, Torquato and 
Jiao [25] found an algorithm which generated jammed 
states over a wide range of packing fractions. (The com- 
monly used Lubachevsky-Stillinger algorithm [26] mainly 
produces jammed states with <pj « 0.635. When the 
number of particles is small, the value of (j)j is reduced 
and the distribution of the range of broadens [27].) We 
attribute the slow dynamics in the system to the squeez- 
ing through narrow necks in the configuration space 
which is required to escape from the vicinity of the initial 
conditions. 

The key to the derivation of our main result, ip = is 
the following picture of the allowed configuration space 
of N hard spheres or disks in a box. Consider spheres. 
As the packing fraction increases beyond <fi ~ 0.59, there 
is marked slowing down of the dynamics of the system; 
r a grows rapidly as the packing fraction increases past 
this density: the spheres are caged and can only escape 
from the cages on the alpha relaxation time scale. We are 
imagining that typical initial conditions place the system 
close to one of the jammed states but that there is a way 
of escaping from this state through a neck in configura- 
tion space. Suppose that the packing fraction <j> is close 
to but below ^Neck- Let 



a = In 



1 



1 - <p/(pNeck 

then from Eq. (7), provided N is large 
r ~ ro exp [aN] . 



/d, a > 0, 



(8) 



(9) 



Note that the form of the equation (i.e. its dependence 
on N) would be hardly altered if a finite fraction of the 
N particles were rattlers. 

The next step in the argument is to apply Eq. (9) in the 
context of the mosaic picture of dynamic heterogeneity 
[28]. The scale of the tiles within the mosaic is £. It 
is pictured that the dynamics of glasses proceeds by a 



series of rearrangements in each tile. Once in one tile the 
atoms have reached a different state; this will sometimes 
trigger state changes in neigboring tiles. Such avalanches 
of movement are commonly observed [28]. To rearrange 
the spheres in a region of size £ a neck will have to be 
passed through in phase space and the time scale for 
doing this will be as given in Eq. (9) but with N replaced 
by the number of spheres in the region of size £ which 
is ~(£/<j) d . Different regions will have different values 
for a, which is in accord with the fact that there is a 
spread of relaxation times associated with the dynamic 
heterogeneities [28]. The alpha relaxation time sets the 
time scale for these processes, so in Eq. (1) it follows 
that if) = d. Note that the mosaic size £ is the size of 
the region on which one has to make re- arrangements in 
order to allow an atom to escape from its surrounding 
cage, so that the timescale r is naturally identified in 
this picture with the alpha relaxation time. 

We should like to thank Chris Fullerton for many useful 
discussions. 
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